Enhanced nonlinear response of Ne 8+ to intense ultrafast x rays 
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We investigate the possible reasons for the discrepancy between the theoretical 
two-photon ionization cross section, ~ 10 -56 cm 4 s, of Ne 8+ obtained within the 
perturbative nonrelativistic framework for monochromatic light [J. Phys. B 34, 
4857 (2001)] and the experimental value, 7 x 10~ 54 cm 4 s, reported in [Phys. Rev. 
Lett. 106, 083002 (2011)] at a photon energy of 1110 eV. To this end, we consider 
Ne 8+ exposed to deterministic and chaotic ensembles of intense x-ray pulses. The 
time-dependent configuration-interaction singles (TDCIS) method is used to quanti- 
tatively describe nonlinear ionization of Ne 8+ induced by coherent intense ultrashort 
x-ray laser pulses. The impact of the bandwidth of a chaotic ensemble of x-ray 
pulses on the effective two-photon ionization cross section is studied within the low- 
est nonvanishing order of perturbation theory. We find that, at a bandwidth of 11 
eV, the effective two-photon ionization cross section of Ne 8+ at a photon energy of 
1110 eV amounts to 5 x 10 -57 and 1.6 x 10~ 55 cm 4 s for a deterministic ensemble 
and a chaotic ensemble, respectively. We show that the enhancement obtained for a 
chaotic ensemble of pulses originates from the presence of the one-photon ls 2 -lsAp 
resonance located at 1127 eV. Using the TDCIS approach, we also show that, for 
currently available radiation intensities, two-photon ionization of a Is electron in 
neutral neon remains less probable than one-photon ionization of a valence electron. 
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I. INTRODUCTION 



Modern highly intense x-ray free-electron lasers (XFELs), such as the Linac Coherent 
Light Source (LCLS) at SLAC National Accelerator Laboratory, Stanford, USA [lj], and the 
SPring-8 Angstrom Compact Free- Electron Laser (SACLA), Harima, Japan ^j, deliver both 
soft and hard x-ray radiation. FLASH at DESY, Hamburg, Germany operates in the 

4|, which is under construction, is 



VUV and soft x-ray regimes, and the European XFEL 



planned to deliver photon energies up to 12 keV. These facilities offer possibilities to explore 
inner-shell electron dynamics and nonlinear response of atoms and molecules to intense x-ray 
radiation (see for example {5- 10|). 

The present theoretical work is triggered by a recent experiment on nonlinear ionization of 
neon atoms performed at the LCLS [6|. The experiment utilized the capability of the LCLS 
to produce unprecedentedly intense x-ray beams, with up to ~10 12 x-ray photons in a ~100 
fs pulse with a peak intensity of ~ 10 17 W/cm 2 . Within a single pulse the initially neutral 
target absorbed multiple photons yielding a variety of ion species in different electronic 
configurations. At aphoton energy of 1110 eV, which is below the i^-shell threshold of 
Ne 8+ , Doumy et al. 6| observed production of hydrogen-like neon, Ne 9+ . The Ne 9+ /Ne 8+ 
ratio was observed to depend quadratically on the peak intensity, which is consistent with 
nonlinear two-photon ionization of Ne 8+ . Nevertheless, the two-photon ionization cross 
section, deduced from this experimental observation with the help of a rate-equation model, 
is 7 x 10~ 54 cm 4 s, which is two orders of magnitude higher than the value obtained within 



perturbation theory 
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In the present paper, we focus on the fo 
observed [6| and theoretically predicted 



lowing points: (i) the discrepancy between the 
111 ] two-photon ionization cross section values of 
Ne 8+ ; (ii) the possibility of two-photon ionization of a Is electron in neutral neon below the 
-fT-shell threshold of neon. 

To describe the nonlinear interaction of Ne 8+ and neutral neon with an intense coher- 
ent ultrashort x-ray pulse, we adopt the time-dependent configuration-interaction singles 
(TDCIS) method — a nonperturbative ab initio multichannel approach 13Ml6j. TDCIS al- 
lows for pulses of arbitrary shape and peak intensity, and provides an intuitive picture of the 
electron dynamics induced by a light pulse of finite duration. Correlation effects between 
the ejected photoelectron and the remaining ion are included via exact treatment of the 
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Coulomb interaction 
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15[. Going beyond the standard single-active electron approximation 



171 ]. the TDCIS model accounts for the coupling between different excitation (ionization) 



channels. 



In our study, we employ the TDCIS method implemented in the xcid code 181 ] . To 



eliminate spurious reflections, which appear when the electronic wave packet reaches the 



boundary of the numerical grid, we app 



complex absorbing potential (CAP) 



19 



y absorbing boundaries through the inclusion of a 



20|. Implemented within the framework of TDCIS, 



the CAP provides a measure for the ionization probability for the outgoing electron. The 
ionization probability, given by the diagonal components of the reduced ion density ma- 
trix (IDM), is used in this work for calculating the generalized two-photon ionization cross 
section. 

For nonlinear lig ht-matter interaction the spectral and temporal shape of the pulse is 



a crucial factor 
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22J. The rate of simultaneous absorption of two photons depends on 



the statistics of the exciting field 23|-|25|. Present XFELs have a coherence time that is 



much shorter than the pulse duration and can be considered as chaotic [26|, [27J]. For a 



chaotic ensemble of pulses 



28|, 



with a finite bandwidth and a short coherence time 
within the lowest nonvanishing order of perturbation theory (LOPT), the effective two- 
photon ionization cross section can be written as a convolution of the monochromatic two- 
photon cross section and the spectral distribution function. We investigate the effect of 
finite coherence time on the two-photon ionization cross section by using a Gaussian spectral 
distribution function. The monochromatic two-photon ionization cross section we calculate 
within the Hartree-Fock-Slater (HFS) model 29], implemented within the xatom code 
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The results to be presented here indicate that the treatment of XFEL radiation as a 
chaotic finite-bandwidth ensemble of pulses, rather than a deterministic ensemble of pulses, 
is likely to be capable of explaining the enhanced two-photon ionization cross section reported 
inRef. Q. 

The paper is organized as follows: we outline the theoretical approaches in Sec. [Ill present 
details on the numerical implementation and the obtained results in Sec. II III and draw 
conclusions in Sec. [TV] Atomic units are used throughout, unless otherwise noted. 
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II. THEORY 

A. Two-photon ionization cross section for a coherent finite pulse 

A detailed description of our implementation of the TDCIS method can be found in 
Ref. {l4j . Briefly, we construct the electronic wave packet in an atom as a linear combination 
of the Hartree-Fock ground state |$ ) an d one-particle-one-hole (lp-lh) excitations |$f), 

|^t)=a (t)|$o) + EE«-(t)|$-), (1) 

i a 

where 

Here, z, j, . . . label orbitals occupied in \$o), whereas unoccupied (virtual) orbitals are 
marked by a, b, . . . The operators & and c pr7 create and annihilate, respectively, electrons 
in a spin orbital of the modified Fock operator -Fcap — F — irjW, which consists of the Fock 
operator F and the CAP in the form —irjW. The spin states are designated with a. In the 
electric dipole approximation, the nonrelativistic Hamiltonian of the atom interacting with 
the x-ray field is given by 

H = Fcap + V c - %f - E HF - S(t)z, (3) 

where Vc stands for the electron-electron Coulomb interaction, Vhf and Fhf are the Hartree- 
Fock mean-field potential and ground-state energy, respectively, z is the dipole operator, and 
S(t) is the electric field of the intense ultrashort laser pulse, which is assumed to be linearly 
polarized along the z axis. By substituting the wave function given by Eq. (TjQ) into the time- 
dependent Schrodinger equation, one gets a set of coupled ordinary differential equations for 
the coefficients oto(t) and 

Using the state \^,t), we construct the reduced density matrix of the residual ion pro- 
duced in the photoionization process, 

p(t) = Tr a [\* 1 t)(*,t\] 1 (4) 

a,b 

where o a b stands for the overlap between eigenfunctions of -F C ap- The CAP is only active 
at large distances from the atom, and, hence, affects only virtual orbitals. Application of 
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the CAP is equivalent to attenuation of the wave packet when it reaches the boundary of 



the numerical grid 32]. Because of the CAP, the norm of the wave packet from Eq. (CQ) as 
well as the norm of the reduced ion density matrix (j2J), are not conserved and decrease as 
ionization proceeds. In order to compensate for this loss of norm in the IDM, one has to 
introduce a correction [^J : 

6 Pij (t) = 2nS ei - e * I dt'^baatit')^')}^-^-^, (6) 

with the £j being the orbital energies and Wb a the matrix elements of the CAP operator 
W. In the limit t — > oo, i.e., after the ionizing pulse is over, a diagonal component of the 
corrected IDM, pi + 5pi(= p„ + 5pu), can be thought of as the excitation probability from 
an occupied orbital %. Under the conditions considered here, the uncorrected p\ vanishes for 
sufficiently long time after the pulse is over, indicating that the photoelectron is completely 
absorbed by the CAP. Conversely, the IDM correction, 5pi, approaches a constant value at 
t — > oo and can be interpreted as the ionization probability of an electron from orbital i. 

The ionization probability per unit time due to direct absorption of iV photons (in s _1 ) 
is given by a^ N \j N where J is the photon flux in number of photons per cm 2 per second. 
This allows for a definition of an effective two-photon ionization cross section for a coherent 
pulse centered at a mean photon energy oo- in with a bandwidth of Au> p 

hm 6pi(t) 

^coh^in, Auj p) = • ( 7 ) 

J dt J(t) 2 

—oo 

The quantities u- m and Acjp enter the right-hand side of Eq. (J7J) implicitly through the IDM 
correction Spi, obtained using the Hamiltonian from Eq. ()3]), and the flux J(t). The definition 
of Eq. (J7|) is valid provided the ground state is not depleted, i.e., in the perturbative limit. 



B. Two-photon ionization cross section for chaotic fields 

When defining the cross section in Eq. (J2J) we assume that the x-ray pulse is well de- 
fined (deterministic). In general, the radiation produced by an XFEL operating in the 
self-amplified spontaneous emission (SASE) regime is chaotic with respect to fluctuations 
in the electric field. The simplest way to account for the XFEL chaoticity is to recall that 
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the iV-photon ionization rate, within the lowest nonvanishing order of perturbation theory 
(LOPT), is proportional to N\J N , which amounts to effective doubling (2!) of the cross sec- 



tion value for two-photon ionization 21] . This factor of two cannot explain the discrepancy 
found in Ref. The most rigorous and accurate way to simulate the experimental situa- 



34. 



35J for the radiation and 



tion would be by introducing an appropriate stochastic model 
solving the TDCIS equations many times using an ensemble of realistic pulses. Afterwards 
one would have to average the results over all members of the ensemble. However, this 
approach is computationally very costly. 

Here, we follow the result of Mollow [25 1 who showed that within the second-order per- 
turbation theory for a field consisting of finite chaotic pulses, the transition rate due to 
two-photon absorption during the pulse can be expressed in terms of the spectral first-order 
field correlation function. In case of a finite-bandwidth field and near an intermediate res- 
onance of the target atom the two-photon ionization cross section for an incoherent pulse 
can be cast in the form: 



die a 



of the LOPT for monochromatic radiation 



^ PT (u)F(uj,uj in ,Auj p ) } (8) 

(2) 



where F(u, u- m , Au p ) is the normalized spectral distribution function and ct^opt ^ s ^ ne resu h 
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4optM = vr(4vrau;) 2 ^5(w / 

/ 



2uj) 



(f\z\l)(l\z\g) 



(9) 



u g + cu - lui + iTi/2 

with a being the fine-structure constant. In Eq. fl9]), |/), \l) and \g) stand for final, in- 
termediate and ground states, respectively. Ti accounts for the natural linewidth of the 
intermediate states ]/); u g and ui denote energies of the ground and intermediate states, re- 
spectively. Note that the factor of 2 in Eq. (p) accounts for the enhancement of two-photon 
absorption from a single-mode chaotic field [211 ] . 

The spectral distribution of a single XFEL pulse is very spiky and random [ 
Averaged over many shots the spectral distribution can be taken as a normalized Gaussian 
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39(], 



F(u,u in , Acj f 



2VhT2 



cxp 



41n2 




(10) 



The result given by Eq. OH]) can be understood as a nonlinear atomic response to a spectral 
range of uncorrelated modes. Here, the atomic response to the individual frequencies is 
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summed incoherently. In contrast, Eq. ([7]) represents nonlinear atomic response to a coherent 
pulse. In the next section, we apply Eqs. (J7J) and (jSJ) to calculate effective two-photon 
ionization cross sections of Ne 8+ in the photon-energy range below its K-edge. 



III. RESULTS AND DISCUSSION 



We start our numerical study with the nonlinear atomic response of Ne 8+ to a determin- 



istic coherent pulse using TDCIS implemented in the xcid code [18|. We obtain converged 
results by using a nonuniform radial grid extending from r=0 to r=80 a.u. with 1000 grid 
points and a pseudospectral-grid parameter £ = 0.461 [ijj]. Under these conditions, there is 
an almost uniform orbital energy spacing of about 0.3 a.u. across a wide energy range (up to 
150 a.u.) for the final states of the outgoing electron. The CAP starts at r=50 a.u. We use 
a CAP strength rj = 0.002 a.u., which makes the energy levels broad enough to describe the 
quasicontinuum. In this range of r], we satisfy the stationarity condition with respect to 77: 
<9[lim 8pi s (t)]/drj = 0, where p\ s denotes the diagonal component of the IDM corresponding 

t— too 

to the Is orbital. The positions of ls 2 -lsnp resonances are obtained with an accuracy of 
0.03 a.u. and the one-photon ionization potential of Ne 8+ equals 43.9 a.u. (1194.1 eV). For 



the comparison, the experimental value of the ionization potential of Ne 8+ is 1195 eV 40]. 
We account for angular momenta of the outgoing electron up to / max = 2. The laser pulse is 
given by S(t) = Eq exp{— 2 In 2(t/r p ) 2 } coa(umt), where r p is the full-width-at-half-maximum 
(FWHM) duration of the pulse intensity 4l|, £ is the peak electric field. 

In Fig. [TJ we show how the Is ionization probability depends on intensity at two different 
photon energies used in the experiment 6j, below (1110 eV) and above (1225 eV) the one- 
photon ionization threshold. For the calculation we use a coherent pulse with a FWHM 
bandwidth of 6 eV. We can see that in double logarithmic scale the slope of the curve 
corresponding to 1110 eV is 2, while that for 1225 eV the slope below saturation is 1. This 
reflects the fact that at 1225 eV, Is ionization is a one-photon process, whereas at 1110 eV, 
it is a two-photon process. Above ~ 3 x 10 18 W/cm 2 , depletion of the ground state becomes 
substantial. 



Doumy et al. 



6| measured the mean photon energy with an uncertainty of several tenths 



of eV and the pulse spectral width was 10 ± 1 eV 42J. In Fig. [21 we show the two-photon 



ionization cross section, calculated using Eq. (J7J) for several pulse durations corresponding 
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FIG. 1: (Color online) Intensity dependence of the ionization probability of Ne 8+ , given by the 
diagonal IDM correction 5p\ s , at photon energies of 1225 eV (above the single-photon ionization 
threshold) and 1110 eV (below the single-photon ionization threshold). A deterministic pulse of 
6-eV bandwidth (FWHM) is used. 

to the FWHM bandwidths of 20, 15, 11, 8 and 6 eV. The peak electric field £ =0.03 a.u. 
was used. Also shown is the cross section 0lq PT (o;) given by Eq. ([9]). For the latter, we use 



the Hartree-Fock-Slater (HFS) model [29|, implemented within the xatom code [30l . l31 |. 
The HFS model positions the intermediate resonances at lower energies than those obtained 
in TDCIS, therefore we shifted the curve 0lq PT (u/) such that the ls 2 -ls4p resonance is at 
the right position of 1127.1 eV. Doumy et al. |6( noticed that in a similar perturbative 
calculation ll| the authors did not account for the ls 2 -ls4p resonance. We have included 
this resonance in both TDCIS and LOPT calculations. However, as we see from Fig. [2l 
neither the inclusion of this resonance nor the finite bandwidth of the radiation pulse taken 
into account in TDCIS can explain the discrepancy of several orders of magnitude between 
the theoretical and experimental values. 

Now, we use Eq. (jSJ) to convolve the monochromatic two-photon ionization cross section 
obtained with Eq. (j^J), with the spectral distribution function given by Eq. (fTUlh and show 
the results in Fig. [3](a). One can see that within the bandwidth, off from the resonances, the 
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FIG. 2: (Color online) Effective two-photon ionization cross section for Ne 8+ . The TDCIS results 
are given by Eq. ([7J) for several different pulse bandwidths (FWHM). The LOPT result is obtained 
using Eq. ([9|). The point at 1110 eV corresponds to the experimental value of 7 x 10~ 54 cm s 
reported in Rcf. 



cross section is substantially enhanced, because the main contribution to the convolution 
in Eq. (JH]) comes from the resonance peaks. Indeed, for a bandwidth of 11 eV the cross 
section at 1110 eV is 1.6 x 10~ 55 cm 4 s, thus is enhanced by at least one and a half orders of 
magnitude with respect to the perturbative result (4x 10~ 57 cm 4 s). In Fig.[3](b), we show the 
relation between the pulse bandwidth Au p and mean photon energy u- m , which is needed for 
the calculated two-photon ionization cross section to reach the experimentally found value 
of 7xl0 -54 cm 4 s. For a bandwidth of 17 eV, the calculated cross section increases up to 
this value at the photon energy of 1110 eV used in the experiment. Thus, our findings 
suggest that the main reasons for the enhanced two-photon ionization cross section of Ne 8+ 
at 1110 eV originate from the proximity of the ls 2 -ls4p resonance, the chaoticity of the 
LCLS radiation, and the finite bandwidth of its pulses. 

In connection with the study of two-photon ionization of core electrons, it is worth men- 
tioning another recent experiment of Young et al. 5|, where direct multiphoton ionization 
of neon was completely shadowed by a sequence of one-photon ionization events. One of the 
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FIG. 3: (Color online) (a) Two-photon ionization cross section for Ne 8+ , given by Eq. ([5]). The 
perturbative result Olqpt °^ © (dotted line) is taken as a reference signal for averaging over 
different bandwidths (FWHM) of the pulses. The point at 1110 eV corresponds to the experimental 



value of 7 x 10 54 cm 4 s reported in Ref. 



6j. (b) Relation between the bandwidth Acj p and the mean 



photon energy ui- m for which the two-photon ionization cross section (j.^^ is fixed at 7x 10 ^ 4 cm 4 s. 



measurements has been done at the photon energy of 800 eV, just below the K-edge, 870 
eV, of neutral neon. In this case, one x-ray photon carries enough energy to ionize valence 
electrons, and therefore the valence-shell electrons are stripped in a sequence of one-photon 
absorption processes. Creation of a ls-shell vacancy is possible only through the absorption 
of two photons. No evidence for this process was detected. 

Using the TDCIS model, we study the possibility of creating a core hole in Ne via si- 
multaneous absorption of two 800-eV photons. The converged result is obtained by using 
a maximum radius of 90 a.u. with 1000 grid points and ( = 0.461. A CAP of strength 
rj = 0.002 starts at 60 a.u. Accounting for angular momenta of the ionized electron up to 
I max = 3 is sufficient. In Fig. |U we show the ionization probabilities of valence and core 
electrons for neutral neon as a function of peak intensity. One can see that at the intensity of 
3 x 10 17 W/cm 2 the probability of ejecting a Is electron is more than 10 2 times smaller than 
that of ejection of a valence electron. With increasing intensity, the relative probability of 
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FIG. 4: (Color online) Intensity dependence of the ionization probability of neutral neon, given 
by the IDM corrections 5p± s for Is electrons and 5p2 S + Ylm ^Pip™. f° r valence electrons, at a 
photon energy of 800 eV (below the one-photon ionization threshold for the ET-shell, but above the 
one-photon ionization threshold for the valence shells). A deterministic pulse of 6-eV bandwidth 
(FWHM) is used. 

Is ionization with respect to valence ionization grows. Nevertheless, this calculation shows 
that direct two-photon processes with ejection of an inner-shell electron never dominate the 
one-photon ionization of valence electrons, even for a pulse as short as r p = 300 as (cor- 
responding to a bandwidth of 6 eV). We confirm the observation of Young et al. |q] that 
multiphoton processes involving inner shell electrons are overshadowed by valence ionization 
as long as the valence electrons are not stripped away. 

IV. CONCLUSION 

In conclusion, we investigated the two-photon ionization cross section of Ne 8+ in the 
vicinity of the ls 2 -lsAp resonance. We presented a strategy for calculating the two-photon 
ionization cross section within the TDCIS framework. However, the TDCIS model, which 
allows for a perfectly coherent radiation pulse, does not explain the enhanced two-photon 



12 

ionization cross section, obtained by Doumy et al. js] at 1110 eV, in spite of the inclusion 



of the ls 2 -ls4p resonance which was missing in Ref. 







The inclusion of the ls 2 -ls4p 



resonance within the LOPT approach for monochromatic light does not explain the experi- 
mental result either. Chaoticity and short coherence time of the XFEL radiation, taken into 
accounted through the spectral distribution function in the cross-section expression obtained 
within LOPT, partially explain the observed enhancement. For the bandwidth of 10 ± 1 eV, 
estimated in the experiment, we obtained an increase of the effective two-photon cross sec- 
tion by a factor of 40 with respect to the perturbative result for monochromatic radiation. 
To explain the experimentally observed value of 7 x 10~ 54 cm 4 s within this framework one 
would need a broader spectral bandwidth (~ 17 eV) or a mean photon energy tuned closer 
to the ls 2 -ls4p resonance. It is also worth noting that some indirect pathways that avoid 
production and two-photon ionization of ground-state Ne 8+ have not been included in the 
rate-equation model used in Ref. 6j. This might have caused the experimental to be 
overestimated. Nevertheless, we believe the ls 2 -ls4p resonance is the key to explain the 
enhanced two-photon ionization cross section of Ne 8+ at 1110 eV, but its influence depends 
strongly on the XFEL spectral density and uncertainties in its mean photon energy. 

From the study of neutral neon performed within the TDCIS framework, we also infer 
that, when available, valence electron stripping due to one-photon ionization dominates 
over two-photon ionization of inner-shell electrons even at intensities far beyond current 
experimental possibilities. 
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